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Abstract 

The effects of the chemical functionalization of a 
carbon nanotube embedded in a 
nanotube/polyethylene composite on the bulk elastic 
properties are presented. Constitutive equations are 
established for both functionalized and non- 
functionalized nanotube composites systems by using 
an equivalent-continuum modeling technique. The 
elastic properties of both composites systems are 
predicted for various nanotube lengths, volume 
fractions, and orientations. The results indicate that 
for the specific composite material considered in this 
study, most of the elastic stiffness constants of the 
functionalized composite are either less than or equal 
to those of the non-functionalized composite. 

Introduction 

Carbon nanotube-reinforced polymer composites 
have the potential to provide order-of-magnitude 
increases in strength and stiffness relative to typical 
carbon- fiber reinforced polymeric composites . 1 The 
remarkable property-to-weight ratio of these materials 
makes them ideal candidates for aerospace-related 
structural applications. The mechanical behavior of 
nanotube/polymer composites depends not only on the 
individual properties of the polymer and the 
nanotubes, but also on the nanotube/polymer 
interaction. 

This polymer/nanotube interaction is assumed to 
be confined to a cylindrical region surrounding the 
nanotube and extending radially away from the 
exterior of the tube. The strength of this 

nanotube/polymer interaction is defined by the ability 
to transfer mechanical load between the nanotube and 
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polymer and has been a matter of significant discussion 
in the literature . 2 " 14 In general, two basic approaches 
have been proposed to improve the strength of the 
interaction. The first approach involves forming a 
strong, non-bonded interaction with the nanotube 
without modifying the nanotube structure. This 
approach assumes a chemical strengthening of the non- 
bonded interactions with the nanotube 15 and the 
improvement of the mechanical connection with the 
nanotube; for example, by wrapping a large polymer 
molecule around the nanotube . 16 ' 17 The second 
approach requires the formation of a chemical covalent 
bond between the nanotube and polymer directly, also 
known as functionalization. There is reasonable 
evidence for the presence of such nanotube/polymer 
covalent bonds . 18 ' 19 Some studies have indicated that 
functionalization can occur through chemical bonds 
added to the nanotube sidewalls . 5 " 8 Addition of small 
organic groups to nanotube sidewalls has been reported 
via in situ - generated diazonium compounds 4 and by 
fluorination of nanotubes 5 ’ 6 followed by alkylation . 7 
Functionalization can also occur through oxidizing the 
nanotube sample to induce the formation hydroxyl or 
carboxylic acid groups at surface defects on the 
nanotube . 9 " 14 This type of functionalization can occur 
with co-polymers , 9 " 11 proteins , 12 organosilanes , 13 and 
metal catalysts . 14 It is known that covalent bonding 
affects the properties of the nanotube itself because the 
formation of a chemical bond with carbon atoms in a 
nanotube interrupts the sp 2 hybridization of the 
nanotube, thereby forming a site that is closer to an sp 3 
hybridized carbon . 20 21 However, analysis and tests 
have not provided a description of how 

functionalization affects the overall mechanical 
properties of the composite. A recent simulation study 
by one of the authors has predicted that for carbon 
nanotube/polyethylene composites there is at least a one 
order-of-magnitude increase in the strength for 
composites with covalent bonding between the 
nanotube and adjacent polymer molecules, relative to 
systems without the covalent bonds . 22 However, the 
hybridization change may weaken the chemical bonds 
in the vicinity of the functionalization. The effect 
should be manifested in the material elastic 
constants . 15 ' 16 It is proposed therefore, that hierarchal 
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analysis and subsequent predictions of macroscopic, 
elastic behavior may provide additional insights into 
the relative merits of nanotube functionalization and 
the relationship of functionalization to 
nanotube/polymer load transfer and strength. 

In the present paper, constitutive models are 
developed for nanotube/polyethylene composites, with 
and without functionalization, to predict the effect of 
functionalization on the elastic mechanical stiffness 
properties. The constitutive models are established 
using a recently developed, hierarchal modeling 
method for predicting the elastic properties of 
nanotube/polymer composites. 23 The constitutive 
models are then used to establish equivalent- 
continuum models of the representative volume 
elements (RVEs) of the functionalized and non- 
functionalized composites. The equivalent-continuum 
model is then combined with a micromechanical 
analysis to predict Young’s modulus and the shear 
modulus, for various nanotube lengths, volume 
fractions, and orientations. Results from these 
analyses are quantified and the predicted properties of 
the functionalized and non-functionalized composite 
systems are compared. 

Equivalent-Continuum Modeling 

The constitutive model of the composite material, 
required by the equivalent-continuum method, uses an 
equilibrium molecular structure obtained with a 
molecular dynamics (MD) simulation. 23 Using a 
hierarchical modeling scheme, the equivalent- 
continuum modeling technique is then used to predict 
the bulk mechanical behavior of nano- structured 
materials. 24 In summary, the hierarchical model 
consists of three major steps. First, a suitable 
representative volume element (RVE) of the nano- 
structured material is chosen based on the geometry of 
the molecular model. The RVE of a typical nano- 
structured material is on the nanometer length scale, 
therefore, the material of the RVE is not continuous, 
but is an assemblage of many atoms. Interaction of 
these atoms is described in terms of molecular- 
mechanics force constants, which are known for most 
atomic structures. 25 Second, an equivalent-truss model 
of the RVE is developed as an intermediate step to link 
the molecular and equivalent-continuum models. 
Each atom in the molecular model is represented by a 
pin-joint, and each truss element represents an atomic 
bonded or non-bonded interaction. The values of the 
modulus of the truss elements are determined using 
the molecular-mechanics force constants. Therefore, 
the total molecular potential energy of the molecular 
model and the strain energy of the equivalent-truss are 
equal for the same loading conditions. Finally, an 
equivalent-continuum model of the RVE is developed 


in which the total strain energy in the molecular and 
equivalent-continuum models, under identical loading 
conditions, is equal. The effective mechanical 
properties, or the effective geometry, of the equivalent- 
continuum is then determined from equating strain 
energies. Each step of the implementation of this 
method for a functionalized and non-functionalized 
nanotube/polyethylene composite is described below. 
Further details of the modeling method may be found in 
the previously cited references. 

Molecular model 

MD simulations were used to determine the 
equilibrium structures of the functionalized and non- 
functionalized nanotube/polyethylene composites (left 
side of Figs. 1 and 2). For both systems, the starting 
configuration for the MD simulation was a crystalline 
polyethylene matrix containing a (10,10) nanotube of 
radius 6.78 A. The matrix contained 176 chains of 42 
CH 2 units. The nanotube and polymer chains were 
replicated across the periodic boundaries of the 
simulation cell, making them infinitely long. In the 
functionalized composite, two polymer chains were 
attached to six carbons of the nanotube by chemical 
linkages consisting of 2 CFE groups. The details of the 
linkages are shown in Fig. 3. In the MD simulations, 
all the covalent chemical bonds in the composites were 
modeled with the hydrocarbon potential developed by 
Brenner, et al. 26 The non-bonded interactions between 
the polymer chains and between the nanotube and the 
polymer chains were modeled with the Lennard-Jones 
potential by using the parameters used previously. 22 27 
The system was then equilibrated at 300 K via the MD 
simulation. The simulations were carried out using DL- 
POLY®, 28 which is a package of molecular simulation 
subroutines. The systems simulated here are similar to 
those used in prior work. 22 

Because the equivalent-continuum modeling 
technique used in this study requires the interaction of 
the atoms in the molecular model to be described in 
terms of molecular-mechanics force constants, the 
interaction of the atoms were defined in terms of bond 
stretching, bond-angle variation, and Van der Waals 
interactions. 23 The force constants used for both 
models are shown in Tables 1-3, and were taken from 
the AMBER force field. 29 A detailed discussion of the 
parameters used in Tables 1-3 may be found 
elsewhere. 23 ' 29 The force constants were assigned with 
respect to the structure indicated in Fig. 3. The atoms 
labeled as CA, CT, and F1C correspond to aromatic 
carbon atoms (sp 2 orbital), aliphatic carbon atoms (sp 3 
orbital), and hydrogen atoms that are attached to 
aliphatic carbon atoms, respectively. In the non- 
functionalized system, all of the carbon atoms in the 
nanotube were modeled as aromatic carbon atoms. 
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Truss model 

To implement the resultant equivalent-truss 
structures, finite element models were developed by 
using ANSYS® 6 30 (middle portion of Figs. 1 and 2). 
Each element (L1NK8) was a three-dimensional pin- 
jointed truss element with six degrees of freedom 
(three displacement components on each end) that 
represented a single atomic interaction. Each node 
corresponded to an atom in the equilibrium structure 
of the molecular model. The Young’s moduli of the 
truss elements were determined as described 
elsewhere. 23 A total of 36,907 elements and 4751 
nodes were used in the model of the non- 
functionalized nanotube, and a total of 36,481 
elements and 4747 nodes were used for the model of 
the functionalized nanotube. 

Continuum model 

With the equivalent-truss structure in place, the 
continuum models were constructed. The geometries 
of the homogeneous, equivalent-continuum RVEs 
were assumed to be cylindrical, similar to that of the 
molecular and truss models (right side of Figs. 1 and 
2). The mechanical properties of the continuum solid 
cylinders were determined by equating the total strain 
energies of the equivalent-truss and equivalent- 
continuum models under identical loading 
conditions. 24 ' 23 The molecular model of the 
functionalized nanotube/polyethylene composite was 
assumed to exhibit orthotropic symmetry. There are 
nine independent material parameters required to 
determine the entire set of elastic constants for an 
orthotropic material and each of the nine parameters 
may be determined from a single boundary condition 
applied to both equivalent-truss and equivalent- 
continuum models. The non-functionalized molecular 
model had transversely-isotropic symmetry, and 
therefore could be completely described by only five 
independent parameters. Flowever, for consistency 
and without loss of accuracy, nine independent 
parameters were determined with the same nine 
boundary conditions used in the functionalized model. 
Once the mechanical properties of the equivalent- 
continuum RVEs were determined, the RVEs were 
assumed to behave in the composite as effective fibers 
and were used in subsequent micromechanical 
analyses. 

Effective-fiber constitutive model 

The constitutive relationship of an orthotropic 
equivalent-continuum RVE (which is referred to as the 
effective fiber throughout the remainder of the paper) 
is 


ct — O 4 p + O 4 p + O 4 p 

U 11 ^ll c ll T '-'12* > 22 T ^13 & 
CT-,-, = C( 2 

CT = ( P 4- f 4 P + O 4 P 

u 33 ^13°11 ' '^23°22 ~ ^33° 

ct — 2f 4 p 

u 23 ^^44^23 

ct = 2r 4 p 

u 13 z -'^55 & 13 

ct = 2 O 4 p 
u 12 z, '^66 & 12 


( 1 ) 


where cr and s i; are the stress and strain components, 
respectively ( ij = 1,2,3), and C[ t are the elastic 
stiffness components of the effective fiber (denoted by 
superscript f) which are written in the usual contracted 
notation ( k,l = 1,..,6). 

Nine independent elastic properties are required to 
describe the overall elastic behavior of an orthotropic 
material, as shown in Eq. (1). For convenience, the 
nine independent elastic properties that can be used to 
describe the overall behavior of the effective fiber are 
the three elastic axial stiffness components, Cj 4 , , C 22 , 
and C 23 ; the three plane-strain bulk moduli, K { : , K{ 3 , 
and K ( 2 ; and the three elastic shear-stiffness 
components, C{ 4 , C( 5 , and C{ 6 . The three plane-strain 
bulk moduli are defined as 

K( } =C{ 2 +Ci+2C{ 3 

Kn = C{ x + C{ 3 + 2C( 2 (2) 

K n=Cn+C( 2 +2C{ 2 

where the subscripts indicate the plane that is subjected 
to a plane-strain deformation. Equation (2) can be 
easily derived in a manner similar to the derivation of 
the transverse bulk modulus for a transversely-isotropic 
material.’ 3 Once the nine independent elastic properties 
are determined, the elastic interaction stiffness 
components can be calculated from the relations in Eq. 
(2), that is, 


ci=\(zL-cL-c&) 

(3) 

C { 2 = l -{K( 2 -C f n -C f 22 ) 

At this point, the values of the nine elastic parameters 
are unknown. These values are determined by applying 
nine identical sets of boundary conditions to the 
equivalent-truss model and the effective fiber, and by 
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subsequently equating the strain energies by adjusting 
the nine independent elastic properties. The boundary 
conditions used were chosen such that each of the nine 
corresponding strain energies of the effective fiber was 
dependent on only one parameter. 

Boundary conditions 

Nine sets of boundary conditions were chosen to 
uniquely determine each of the nine independent 
elastic properties. The displacements applied at the 
boundaries of the RVE are generalized by 

u,( B ) = ZyXj (4) 

where B is the bounding surface, Xj is defined in Figs. 
1 and 2, and the summation convention associated 
with repeated indices is used. The total strain energy 
of the effective fiber is 


analysis by using the mechanical properties of the 
effective fiber and the bulk polymer matrix material. 
For the composites considered in this study, the 
polymer molecules that were near the 

polymer/nanotube interface were included in the 
effective fiber, and it was assumed that the matrix 
polymer surrounding the effective fiber had mechanical 
properties equal to those of bulk polyethylene. All 
relative movement and interaction of the polymer 
chains with respect to each other was modeled at the 
molecular level. This motion and interaction is 

therefore indirectly considered in the subsequent 
determination of the properties for the effective fibers, 
and it is therefore assumed that perfect bonding exists 
between the effective fiber and the surrounding 
polymer matrix in the micromechanics analysis. 

The micromechanics-based Mori-Tanaka 
method 31 32 was used to predict the elastic mechanical 
properties of the composite material considered herein. 
For this method, the overall elastic-stiffness tensor of 
the composite containing orthotropic effective fibers 
embedded in a matrix material is 


where V is the volume of the effective fiber and 
includes the boundary region surrounding the 
nanotube. The boundary conditions and strain 
energies for each of the nine independent elastic 
properties are listed in Tables 4-6. Unspecified strain 
components in Tables 4-6 are zero valued. 

Boundary region 

The displacements specified in Tables 4-6 were 
applied to each node in the boundary region of the 
equivalent-truss model (indicated in Figs. 1 and 2), 
and the corresponding strain energies were calculated 
by summing the strain energies of each individual 
truss member in the RVE. 

It was assumed that the thickness of the boundary 
region was equal to the maximum distance for which a 
positive-definite relationship exists between the force 
and displacement of all of the atom types. It was also 
assumed that the contribution of the energies 
associated with van der Waals forces between atoms 
with a larger distance than this maximum were 
relatively small and could be neglected. 

Based on the repeating unit used in the MD 
simulation, the effective fiber radius and length were 
chosen to be 1.1 nm 4.3 nm, respectively. The 
calculated values of the nine independent parameters 
for the functionalized and non-functionalized effective 
fibers are listed in Table 7. 

Micromechanical modeling 

Constitutive models of the effective fiber/polymer 
composite were developed with a micromechanical 


C = (c m C m + C/ (C'T' ))(c m I + C/ (T 2 ' )) ' (6) 


where Cf and c m are the effective fiber and matrix 
volume fractions, respectively, Cf and C m are the 
stiffness tensors of the effective fiber and matrix, 
respectively, I is the identity tensor, the angle-brackets 
indicate an effective-fiber orientation average, and T 1 is 
the dilute strain-concentration tensor of the effective 
fibers, and is given by 


T f = 


I + S f (C"‘ ) ' (c f 



( 7 ) 


where Sf is the Eshelby tensor. 33 For unidirectionally- 
aligned effective fibers, the orientation averaging in Eq. 
(6) is not necessary, and the resulting elastic stiffness 
components of the composite have orthotropic 
symmetry. For three-dimensional randomly oriented 
effective fibers, the orientation average of a tensor. A, 
is 


(A) = (4_> = (k-^]( 5,8„) (g> 

+ u(8.8. +5.5 ) 

where ij,m,n = 1,2,3; the indicial summation 

convention is used; 8,y is the Kronecker delta; and 
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K = \ A ‘W 

_i 4 1 

• 10 'C* 3 " "J 

Therefore, from Eqs. (8) and (9), (A) is isotropic. 

For the effective fiber/polymer composite 
considered in the present study, the elastic stiffness 
components, volume fraction, length, and orientation 
of the effective fiber were used for the inclusion 
properties in Eq. (6). The effective fibers were 
assumed to have an spheroidal geometry for the 
Eshelby tensor, so that any length of effective fiber 
could be modeled. The effective-fiber length was 
adjusted by varying the inclusion aspect ratio while 
holding the diameter constant. While the nanotube and 
effective-fiber lengths are equivalent, the nanotube 
volume fraction was determined to be 52.9% of the 
effective-fiber volume fraction. This value was 
calculated by assuming the nanotube volume is 
defined as the total space occupied by the nanotube, 
including half of the unoccupied space between the 
nanotube and polymer. It was assumed that the 
polyethylene matrix was isotropic, with a 
representative Young’s modulus and Poisson’s ratio of 
0.9 GPa and 0.3, respectively. The overall composite 
stiffness was calculated for various effective-fiber 
lengths and volume fractions for both cases of 
perfectly aligned and three-dimensional randomly- 
oriented effective fibers. 

Results 


In this section, the moduli of the nanotube 
composites are presented in terms of nanotube length, 
volume fraction, and orientation. While the shear 
moduli of the composite material (G 44 , G 55 , G 66 ) are 
simply equal to the shear-stiffness components (C 44 , 
C 55 , C 66 ), the Young’s moduli (E 1; E 2 , E 3 ) were 
calculated by using the components of the compliance 
tensor of the composite material, which was 
determined by inverting the composite stiffness tensor, 
C . 34 The subscripts of the Young’s moduli and shear 
moduli indicate the principal direction associated with 
the quantity, similar to Eq. (1). Therefore, Ei is the 
longitudinal Young’s modulus (parallel to the 
nanotube), E 2 and E 3 are transverse Young’s moduli, 
G 44 is the transverse shear modulus, and G 55 and G w , 
are longitudinal shear moduli. Since the random 
composites are isotropic, their mechanical properties 
are completely described by their Young’s modulus, E 
(E = E[ = E 2 = E 3 ), and shear modulus, G (G = G 44 = 
G 55 = G 66 ). 


The Young’s modulus of the nanotube composite is 
shown in Fig. 4 as a function of nanotube length, for a 
1% nanotube volume fraction. The longitudinal 
Young’s modulus of the aligned composite and the 
Young’s modulus of the random composite have the 
expected non-linear dependence on the nanotube 
length. An increase in the Young’s modulus with 
respect to an increase in nanotube length is expected to 
correspond to a relative increase in load transfer 
efficiency between the nanotube and the surrounding 
polymer. The data in Fig. 4 indicate that at a nanotube 
length of about 400 nm, the efficiency of load transfer 
is nearly maximized. Further increases in nanotube 
length result in relatively small increases in Young’s 
modulus for a given volume fraction. As the nanotube 
length becomes greater than approximately 100 nm, the 
Young’s modulus for the non-functionalized composite 
becomes larger than that of the functionalized 
composite for the random composite and longitudinal 
deformation of the aligned composite. By contrast, the 
transverse Young’s modulus of the aligned composite 
has no dependence on nanotube length. Also, there is 
no effect of functionalization on the transverse 
deformation of the aligned composite. 

The shear modulus of the nanotube composite is 
shown in Fig. 5 as a function of nanotube length. The 
shear modulus of the random composite shows an 
expected non-linear dependence on nanotube length, 
whereas the shear modulus of the aligned composite has 
no dependence on nanotube-length. Again, the data 
indicate that the efficiency of load transfer, as 
evidenced by a modulus change, is nearly maximized at 
a nanotube length of about 400 nm. For the random 
composite, the shear modulus of the non-functionalized 
system increases over that of the functionalized system 
as the nanotube length is increased. There is no 
significant effect of functionalization on the shear 
modulus of the aligned composite as the nanotube 
length is increased. As the nanotube length approaches 
the length of the effective fiber, 4.3 nm, the shear 
modulus is significantly higher than the shear modulus 
of the polyethylene matrix, 0.35 GPa (calculated using 
the isotropic polyethylene properties given in the 
previous section), which indicates the large effect of 
small inclusions on the shear modulus at fixed 
concentrations. 

The Young’s moduli of the random and aligned 
composites are shown in Figs. 6 and 7 as a function of 
nanotube volume fraction, for a constant nanotube 
length of 400 nm. At this nanotube length, a change in 
the longitudinal Young’s modulus with respect to 
volume fraction of the aligned composite is nearly 
linear, resembling a linear rule-of-mixtures relationship. 
Over the complete range of volume fraction. Young’s 
modulus for the non-functionalized composite is greater 
than that of the functionalized composite for both the 
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random composite and longitudinal deformation of the 
aligned composite. Conversely, for the transverse 
deformation of the aligned composite, the Young’s 
modulus of the aligned, functionalized composite is 
greater than the non-functionalized material. 

The shear modulus of the nanotube composite is 
shown in Fig. 8 as a function of nanotube volume 
fraction. As the nanotube volume fraction increases, 
the shear modulus for the non-functionalized 
composite becomes greater than that of the 
functionalized composite for the random-nanotube 
orientation. There is no significant effect of the 
functionalization on the shear modulus of the aligned 
composite for the given range of nanotube volume 
fractions. 

Summary and Conclusions 

In this study, the bulk composite elastic properties 
of a nanotube/polyethylene composite have been 
predicted for both functionalized and non-functionalized 
nanotubes. It was assumed that functionalization would 
effect the load transfer between the nanotube and the 
surrounding matrix, and that changes in predicted elastic 
stiffness would be directly related to the efficiency of the 
load-transfer mechanism. The analysis method 
establishes constitutive equations for both functionalized 
and non-functionalized nanotube composites systems by 
using a hierarchical equivalent-continuum modeling 
technique that predicts elastic bulk behavior, using 
intrinsic properties developed through molecular 
dynamics simulations. The elastic properties of the 
functionalized and non-functionalized composites 
systems have been examined for various nanotube 
lengths, volume fractions, and both random and aligned 
nanotube orientations. 

For a fixed nanotube volume fraction of 1% and 
various nanotube lengths, the Young’s modulus of the 
random composite, the shear modulus of the random 
composite, and the longitudinal Young’s modulus of 
the aligned composite have been shown to decrease up 
to 10% when the nanotube is functionalized. For a 
fixed nanotube length of 400 nm and various nanotube 
volume fractions, the longitudinal Young’s modulus of 
the aligned composite, the Young’s modulus of the 
random composite, and the shear modulus of the 
random composite have also been shown to decrease 
up to 12% when the nanotube is functionalized. 
However, under these conditions, the transverse 
Young’s modulus of the aligned composite has shown 
an increase of up to 24% when the nanotube is 
functionalized. 

The primary implication of these results is that 
although chemical functionalization of carbon 
nanotubes has been considered a possible means to 
increase load-transfer efficiency in a 


nanotube/polymer composite, this functionalization has 
in fact degraded the macroscopic elastic stiffness of the 
composite considered in this study. In general, occurs 
for both random composites and aligned composites 
over a wide range of nanotube volume fractions and 
lengths. 
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Table 1. Bond-stretching parameters 


Bond 

stretching 

P (A) 

K p 

(kcal/mole/A 2 ) 

CA-CA 

1.400 

469 

CA-CT 

1.510 

317 

CT-CT 

1.526 

310 

CT-HC 

1.090 

340 


Table 2. Bond-angle variation parameters 


Bond-angle 

variation 

© (deg) 

K° 

(kcal/mole/rad 2 ) 

CA-CA-CA 

120.0 

63 

CA-CA-CT 

120.0 

70 

CA-CT-CA 

120.0 

63 

CA-CT-CT 

114.0 

63 

CT-CT-CT 

109.5 

40 

CT-CT-HC 

109.5 

50 

HC-CT-HC 

109.5 

35 


Table 3. Van der Waals interaction parameters 


Van der Waals 
interaction 

D l 

(kcal/mole) 

P J (A) 

CA 

0.0860 

3.816 

CT 

0.1094 

3.816 

HC 

0.0157 

2.974 
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Table 4. Boundary conditions for axial stiffness 
components 


Property 

Boundary 

conditions 

Boundary 

displacements 

Strain energy 

C f 

Mi 

s n =e 

M| (f?) = ex l 
u 2 (B) = 0 
n 3 ( B ) = 0 

A < = V -Ctf 

c f 

^22 

s 22 — s 

M, (5) = 0 
u 2 (f?) = ex 2 
u 3 (f?) = 0 

a r = v -cu 

C f 

M3 

833 =e 

Mj (5) = 0 

« 2 (5) = 0 
u 2 (f?) = ex 2 



Table 5. Boundary conditions for plane-strain bulk 

moduli 


Property 

Boundary 

conditions 

Boundary 

displacements 

Strain energy 


S 22 = S 33 = e 

Mj (5) = 0 

u 2 (B) = ex 2 
u 2 ( B ) = ex 2 

A f =-K'e 2 
2 


S 11 = S 33 = e 

u x (B) = ex l 
u 2 ( B ) = 0 
u 2 (f?) = ex 2 

A f =-KU 2 

2 13 

A 

S 11 = S 22 = e 

u ] (f?) = ex 1 
u 2 ( B ) = ex 2 
u 2 ( B ) = 0 

A f =-K{ 2 e 2 
2 12 
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Table 6. Boundary conditions for shear-stiffness 
components 


Property 

Boundary 

conditions 

Boundary 

displacements 

Strain energy 

C f 

'“'44 

^23 = Y/2 

M| (5) = 0 
M 2 { B ) = ^h 

M 3 ( B ) = \ X 2 

A'=|c'y ! 

C f 

^55 

£i3 = y/ 2 

M l( 5 ) = ^ X 3 

w 2 (5) = 0 

U 3 ( B ) = ^ X I 

A / =^cir 

c f 

^66 

£i2 = y/ 2 

U A B ) = \ X 2 

U l( B ) = \ X 1 

n 3 (5) = 0 

A^^cir 2 


Table 7. Effective-Fiber Properties (GPa) 



Non- functionalized 
nanotube 

Functionalized 

nanotube 

c f 

'-'11 

548.7 

487.7 

c f 

^22 

16.8 

24.5 

Q 7 3 

16.5 

20.6 

K{. \ 

59.0 

77.9 

k { 3 

597.1 

548.2 

K \2 

597.0 

554.6 

c f 

'“'44 

7.1 

12.7 

c f 

^55 

144.0 

155.4 

C f 

^66 

144.9 

137.0 
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Figure 1. Equivalent-continuum modeling of non-functionalized system 
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Figure 2. Equivalent-continuum modeling of functionalized system 
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Figure 3. Modeled atomic structure of the functionalized nanotube/polyethylene composite 



Figure 4. Young’s moduli of both composite systems vs. nanotube length for a 1% nanotube volume fraction 
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Figure 5. Shear moduli of both composite systems vs. nanotube length for a 1% nanotube volume fraction 
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Figure 6. Longitudinal Young’s moduli of both composite systems vs. nanotube volume fraction for a nanotube 
length of 400 nm 
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Figure 7. Transverse Young’s moduli of both composite systems vs. nanotube volume fraction for a nanotube length 
of 400 nm 



Figure 8. Shear moduli of both composite systems vs. nanotube volume fraction for a nanotube length of 400 nm 
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